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Abstract 



q ! The relativistic six-quark equations including the u, d quarks and antiquarks are found. The 
CN ; nonstrange baryonia BB are constructed without the mixing of the quarks and antiquarks. The 
P-r relativistic six-quark amplitudes of the baryonia are calculated. The poles of these amplitudes 
determine the masses of baryonia. 15 masses of baryonia are predicted. The mass of baryonium 
with the spin-parity J p = 0~ M = 1835 MeV is used as a fit. 
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BES Collaboration observed a significant threshold enhancement of pp mass spectrum in the 
>• 1 radiative decay J/ip — > ^pp [1]. Recently BES Collaboration reported the results on X(1835) in 
^ ■ the J/ij) — > 7 7/7r + 7r~ channel [2]. Under the strong assumption that the pp threshold enhancement 
CD ! and X(1835) are the same resonance, Zhu and Gao suggested Y(1835) could be a J PC = 0~ + 
I I G = + pp baryonium [3]. 

Theoretical investigations of baryon-antibaryon bound states date back to the proposal of 
Fermi and Yang [4] to make the pion out of a nucleon-antinucleon pair. The model of Nambu 
and Jona-Lasinio [5] which is constructed to give a nearly zero-mass pion as a fermion-antifermion 
bound state, also has a scalar resonance of twice the fermion mass. Enhancement in the baryon- 
antibaryon channel near the threshold are expected on the basis of duality arguments [6 - 8] and 
by comparison with the systematic of resonance formation in meson-meson and meson-baryon 
& \ channels [9]. A historical survey of bound states or resonances coupled to the nucleon-antinucleon 
channel is given in Ref. [10]. Gluonic states can couple to baryon-antibaryon channels of ap- 
propriate spin and parity. The discussions of B decays involving baryon-antibaryon pairs include 
Refs. [11-15]. 

Theoretical work speculated many possibilities for the enhancement such as the t-channel pion 
exchange, some kind of threshold kinematical effects, as new resonance below threshold or pp 
bound state [16 - 23]. 

In a series of papers [24 - 28] a method has been developed which is convenient for analysing 
relativistic three-hadron systems. The physics of the three-hadron system can be described by 
means of a pair interaction between the particles. There are three isobar channels, each of which 
consists of a two-particle isobar and the third particle. The presence of the isobar representation 
together with the condition of unitarity in the pair energies and of analyticity leads to a system of 
integral equations in a single variable. Their solution makes it possible to describe the interaction 
of the produced particles in three-hadron systems. 
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In Ref. [29] a representation of the Faddeev equation in the form of a dispersion relation in the 
pair energy of the two interacting particles was used. This was found to be convenient in order 
to obtain an approximate solution of the Faddeev equation by a method based on extraction of 
the leading singularities of the amplitude. With a rather crude approximation of the low-energy 
NN interaction a relatively good description of the form factor of tritium (helium-3) at low q 2 
was obtained. 

In our papers [30 - 32] relativistic generalization of the three-body Faddeev equations was 
obtained in the form of dispersion relations in the pair energy of two interacting quarks. The 
mass spectrum of S-wave baryons including u, d, s quarks was calculated by a method based on 
isolating the leading singularities in the amplitude. We searched for the approximate solution 
of integral three-quark equations by taking into account two-particle and triangle singularities, 
all the weaker ones being neglected. If we considered such approximation, which corresponds to 
taking into account two-body and triangle singularities, and defined all the smooth functions of 
the subenergy variables (as compared with the singular part of the amplitude) in the middle point 
of the physical region of Dalitz-plot, then the problem was reduced to the one of solving a system 
of simple algebraic equations. 

In the previous paper [35] the relativistic six-quark equations are found in the framework of 
coupled-channel formalism. The dynamical mixing between the subamplitudes of hexaquark are 
considered. The six-quark amplitudes of dibaryons are calculated. The poles of these amplitudes 
determine the masses of dibaryons. We calculated the contribution of six-quark subamplitudes to 
the hexaquark amplitudes. 

In the present paper the relativistic six-quark equations including u, d quarks and antiquarks 
are found. The nonstrange barionia BB are constructed without the mixing of the quarks and 
antiquarks. The relativistic six-quark amplitudes of the baryonia are calculated. The poles of 
these amplitudes determine the masses of baryonia. In Sec. II we briefly discuss the relativistic 
Faddeev approach. The relativistic three-quark equations are constructed in the form of the 
dispersion relations over the two-body subenergy. The approximate solution of these equations 
using the method based on the extraction of leading singularities of the amplitude are obtained. 
We calculated the mass spectrum of S-wave baryons with J p = | + , | + (Table I). In Sec. Ill the six- 
quark amplitudes of baryonia are constructed. The dynamical mixing between the subamplitudes 
of baryonia are considered. The relativistic six-quark equations are obtained in the form of the 
dispersion relations over the two-body subenergy. The approximate solutions of these equations 
using the method based on the extraction of leading singularities of the amplitude are obtained. 
Sec. IV is devoted to the calculation results for the baryonia mass spectrum and the contributions 
of subamplitudes to the baryonia amplitude (Tables II, III, IV). In conclusion, the status of the 
considered model is discussed. 

II. Brief introduction of relativistic Faddeev equations. 

We consider the derivation of the relativistic generalization of the Faddeev equation for the 
example of the A-isobar (J p = | + ). This is convenient because the spin-flavour part of the wave 
function of the A-isobar contains only nonstrange quarks and pair interactions with the quantum 
numbers of a J p = 1 + diquark (in the color state 3 C ). The 3q baryon state A is constructed as 
color singlet. Suppose that there is a A-isobar current which produces three u quarks (Fig. la). 
Successive pair interactions lead to the diagrams shown in Fig. lb-lf. These diagrams can be 
grouped according to which of the three quark pairs undergoes the last interaction i.e., the total 
amplitude can be represented as a sum of diagrams. Taking into account the equality of all pair 
interactions of nonstrange quarks in the state with J p = 1 + , we obtain the corresponding equation 
for the amplitudes: 
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A(s, S12, si 3 , S23) = A + Al(s, S12) + Ai(s, si 3 ) + Ai(s, S23) • (1) 

Here, the s,ik are the pair energies of particles 1, 2 and 3, and s is the total energy of the system. 
Using the diagrams of Fig. 1, it is easy to write down a graphical equation for the function 
Ai(s, S12) (Fig. 2). To write down a concrete equations for the function Ai(s, si 2 ) we must specify 
the amplitude of the pair interaction of the quarks. We write the amplitude of the interaction of 
two quarks in the state J p = 1 + in the form: 

/ \ G\{s 12 ) 



'12 

Bl (s 12 )= J ^k P^2)Gi(s' 12 ) 

J 71 S' 12 -S 12 

4m 2 




M^) = (^t^)F 1 — ^1 ■ (D 

Here G\(s 12 ) is the vertex function of a diquark with J p = 1 + . Bi(s 12 ) is the Chew-Mandelstam 
function [34], and pi(si 2 ) is the phase spaces for a diquark with J p = 1 + . 

The pair quarks amplitudes qq — > qq are calculated in the framework of the dispersion N/D 
method with the input four-fermion interaction [35, 36] with the quantum numbers of the gluon 
[37, 38]. 

The four-quark interaction is considered as an input: 

g v (qXI f -f„q) 2 + 2 g$ (qXIfy^q) (sX-y^s) + g { y S) (sX-y^s) 2 . (5) 

Here If is the unity matrix in the flavor space (u,d). X are the color Gell-Mann matrices. Di- 
mensional constants of the four-fermion interaction gy, gy^ and gy S ^ are parameters of the model. 
At gy = gy^ = gy S ^ the flavor SU (3) / symmetry occurs. The strange quark violates the flavor 
SU (3) / symmetry. In order to avoid additional violation parameters we introduce the scale of the 
dimensional parameters [38]: 

m 2 _ (m + m s ) 2 (s) _ m 2 {ss) 
9 ~ -^9v ^5 9v ~ -^9v , (6) 



A 4A(^fc) 
(rrii + m k ) 2 ' 

Here mi and are the quark masses in the intermediate state of the quark loop. Dimensionless 
parameters g and A are supposed to be constants which are independent of the quark interaction 
type. The applicability of Eq. (5) is verified by the success of De Rujula-Georgi-Glashow quark 
model [37], where only the short-range part of Breit potential connected with the gluon exchange 
is responsible for the mass splitting in hadron multiplets. 

In the case under discussion the interacting pairs of particles do not form bound states. Thefore, 
the integration in the dispersion integral (7) run from 4m 2 to 00. The equation corresponding to 
Fig. 2 can be written in the form: 
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^l(s,Sl2) = 



Algl(gl2) gl(glg) 7 <fei2 Pl(Si 2 ) r / / 

l-£i(s 12 ) i vr s' 12 -s 12 11 " J 

4m 2 



12 ; 



— [A 1 {s,s' 13 )+A 1 (s,s' 23 )} 



(7) 



In Eq. (7) z is the cosine of the angle between the relative momentum of particles 1 and 2 in 
the intermediate state and the momentum of the third particle in the final state in the c.m.s. of 
the particles 1 and 2. In our case of equal mass of the quarks 1, 2 and 3, s' 13 and s' 12 are related 
by the equation (8) (See Ref. [39]) 



s 13 = 2m — 



{s' 12 + m 2 - s) z_ 
2 2\ 



<A ' 2 lm %' 12 -(V-s + m)i)(s' 12 -(V~s-m) 2 ). (8) 



s 12 



The expression for s 23 is similar to (8) with the replacement z —> —z. This makes it possible 
to replace [A^s, s' 13 ) + A^s, s' 23 )\ in (7) by 2A 1 (s, s' 13 ). 

From the amplitude Ai(s, Si 2 ) we shall extract the singularities of the diquark amplitude: 



Ai(s, si 2 ) 



ai(s,Si 2 )-B 1 (s 12 ) 



1 - B^sn) 

The equation for the reduced amplitude a±(s, Si 2 ) can be written as 



(9) 



Am 2 -1 



(h 2ai(s,g / 13 ).Bi(s / 1 3) 
2 l-5i(s' 13 ) 



(10) 



The next step is to include into (10) a cutoff at large s' 12 . This cutoff is needed to approximate 
the contribution of the interaction at short distances. In this connection we shall rewrite Eq. (10) 

as 
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Ql(s , Sl2) = a + — ^ — - J ^e ( A - s< 12 )^L Gl f % 2a f . 

Bi{s 12 ) J o vr s' 12 -s 12 J 2 1-Bi(s 13 ) 

In Eq. (11) we have chosen a hard cutoff. However, we can also use a soft cutoff, for instance 
Gi(s[ 2 ) = G\ exp [— <ySl2 ~ J ^ n - ^ , which leaves the results of calculations of the mass spectrum 
essentially unchanged. 

The construction of the approximate solution of Eq. (11) is based on extraction of the leading 
singularities are close to the region Sik ~ 4m 2 . The structure of the singularities of amplitudes 
with a different number of rescattering (Fig. 1) is the following [39]. The strongest singularities 
in arise from pair rescatterings of quarks: square-root singularity corresponding to a threshold 
and pole singularities corresponding to bound states (on the first sheet in the case of real bound 
states, and on the second sheet in the case of virtual bound states). The diagrams of Figs, lb and 
lc have only these two-particle singularities. In addition to two-particle singularities diagrams of 
Figs. Id and le have their own specific triangle singularities. The diagram of Figs. If describes a 
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larger number of three-particle singularities. In addition to singularities of triangle type it contains 
other weaker singularities. Such a classification of singularities makes it possible to search for an 
approximate solution of Eq. (11), taking into account a definite number of leading singularities 
and neglecting the weaker ones. We use the approximation in which the singularity corresponding 
to a single interaction of all three particles, the triangle singularity, is taken into account. 

For fixed values of s and s' 12 the integration is carried out over the region of the variable s' 13 
corresponding to a physical transition of the current into three quarks (the physical region of 
Dalitz plot). It is convenient to take the central point of this region, corresponding to z — 0, to 
determinate the function ai(s, su) and also the Chew-Mandelstam function -Bi(si 2 ) at the point 
Si 2 = s = § + m 2 . Then the equation for the A isobar takes the form: 

&i(s, s ) = A + h,i(s, s ) • 2 «i(s, s ) , (12) 



/.,.(..*)= J d -^J^L Gl J* — * (13) 

J it s' 12 -s 12 J 2 1-Bi(s' 13 ) 

4m 2 —1 

We can obtain an approximate solution of Eq. (14) 

a 1 (s,s ) = A[l- 2/^(5, so)]" 1 . (14) 

The function ii,i(s, s ) takes into account correctly the singularities corresponding to the fact 
that all propagators of triangle diagrams like those of Figs. Id and le reduce to zero. The right- 
hand side of (14) may have a pole in s, which corresponds to a bound state of the three quarks. 
The choice of the cutoff A makes it possible to fix the value of the mass of the A isobar. 

Baryons of ^-wave multiplets have a completely symmetric spin-flavor part of the wave function, 
and spin | corresponds to the decuplet which has a symmetric flavor part of the wave function. 
Octet states have spin | and a mixed symmetry of the flavor function. 

In analogy with the case of the A isobar we can obtain the rescattering amplitudes for all S- 
wave states with J p = | + , which include quarks of various flavors. These amplitudes will satisfy 
systems of integral equations. In considering the J p = | + octet we must include the integration 
of the quarks in the + and 1 + states (in the colour state 3 C ). Including all possible rescattering of 
each pair of quarks and grouping the terms according to the final states of the particles, we obtain 
the amplitudes Aq and A\, which satisfy the corresponding systems of integral equations. If we 
choose the approximation in which two-particle and triangle singularities are taken into account, 
and if all functions which depend on the physical region of the Dalitz plot, the problem of solving 
the system of integral equations reduces to one of solving simple algebraic equations. 

In our calculation the quark masses m u = = m and m s are not uniquely determined. 
In order to fix m and m s anyhow, we make the simple assumption that m = |?tia(1232) m = 
|mn(1672). The strange quark breaks the flavor 577(3)/ symmetry (6). 

In Ref. [32] we consider two versions of calculations. If the first version the S77(3) / symmetry is 
broken by the scale shift of the dimensional parameters. A single cutoff parameter in pair energy 
is introduced for all diquark states Ai = 12.2. 

In the Table I the calculated masses of the S'-wave baryons are shown. In the first version we use 
only three parameters: the subenergy cutoff A and the vertex function g , g±, which corresponds 
to the quark-quark interaction in + and 1 + states. In this case the mass values of strange 
baryons with J p = | + are less than the experimental ones. This means that the contribution 
color-magnetic is too large. In the second version we introduce four parameters: cutoff A , Ai and 
the vertex function g , g x . We decrease the color- magnetic interaction in + strange channels and 
calculated mass values of two baryonic multiplets J p = | + , | + are in good agreement with the 
experimental data [40]. 
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The essential difference between E and A is the spin of the lighter diquark. The model explains 
both the sign and magnitude of this mass splitting. 

The suggested method of the approximate solution of the relativistic three-quark equations 
allows us to calculate the S'-wave baryons mass spectrum. The interaction constants, determined 
the baryons spectrum in our model, are similar to ones in the bootstrap quark model of S'-wave 
mesons [38]. The diquark interaction forces are defined by the gluon exchange. The relative 
contribution of the instanton-induced interaction is less than that with the gluon exchange. This 
is the consequence of l/iV c -expansion [38]. 

The gluon exchange corresponds to the color-magnetic interaction, which is responsible for the 
spin-spin splitting in the hadron models. The sign of the color-magnetic term is such as to made 
any baryon of spin | heavier than its spin-| counterpart (containing the same flavors). 

We manage with quarks as with real particles. However, in the soft region, the quark diagrams 
should be treated as spectral integrals over quark masses with the spectral density p(m 2 ): the 
integration over quark masses in the amplitudes puts away the quark singularities and introduced 
the hadron ones. One can believe that the approximation: 

p(m 2 ) ->■ 5(m 2 - m 2 q ) (15) 

could be possible for the low-lying hadrons (here m q is the "mass" of the constituent quark). 

We hope the approach given by (15) is sufficiently good for the calculation of the low-lying 
baryons being carried out here. The problem of distribution over quark masses is important when 
one considers that the high-excited states need spectral studies. 

III. Six-quark amplitudes of the baryonia. 

We derive the relativistic six-quark equations in the framework of the dispersion relation tech- 
nique. We use only planar diagrams; the other diagrams due to the rules of 1/N C expansion [41 
- 43] are neglected. The current generates a six-quark system. The correct equations for the 
amplitude are obtained by taking into account all possible subamplitudes. It corresponds to the 
division of complete system into subsystems with a smaller number of particles. Then one should 
represent a six-particle amplitude as a sum of 15 subamplitudes: 

^=EV (is) 

i,j=l 

This defines the division of the diagrams into groups according to the certain pair interaction 
of particles. The total amplitude can be represented graphically as a sum of diagrams. We need 
to consider only one group of diagrams and the amplitude corresponding to them, for example 
Ai2- We shall consider the derivation of the relativistic generalization of the Faddeev-Yakubovsky 
approach. 

In our case the low-lying baryonia are considered. We take into account the pairwise interaction 
of all quarks and antiquarks in the baryonia. 

For instance, we consider the l uu -diquarks with spin-parity J p = 1 + for the baryonium con- 
tent uuuddd (Fig. 3). The set of diagrams associated with the amplitude A 12 can further 
be broken down into five groups corresponding to subamplitudes: A\ uu (s, S12345, S1234, S123, S12), 
A\ {s, S12345, S1234, S123, S12J, A[ (s, S12345, S1234, S123, S12J, baryonium A L 2 {s, S12345, S1234, S12, S34 
A3"" 1 " 1 (s, S12345, S12, S34, S56). Here is the two-particle subenergy squared, Sijk corresponds 
to the energy squared of particles i, j, k, s^i is the energy squared of particles i, j, k, I, Sijkim 
corresponds to the energy squared of particles i, j, k, I, m and s is the system total energy squared. 

The amplitude A\ uu (s, S12345, S1234, S123, s i2 ) consists of the five color sub-structures: the diquark 
yiu - m co \ OT s tate 3 C , the quark u in the color state 3 C , the three d antiquarks each in color 
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state 3 C , therefore we obtain 3 C x 3 C = l c + 8 C , 3 C x 3 C x 3 C = l c + 8 C + 8 C + 10*. Then we consider 
the total color singlet. 

The baryonium amplitude Alf™ 1 (s, S12345, S1234, S12, S34) contains the diquark in the color state 
3 C , antidiquark l dd in the color state 3 C , w-quark in the color state 3 C , <i-antiquark in color state 
3 C . We use the following equations: 3 C x 3 C = l c + 8 C , 3 C x 3 C = l c + 8 C . Then the baryonium 
amplitude is the total singlet. 

The amplitude A\ UU1U 1 (s, S12345, S12, S34, s$q) consists of the diquark l uu in the color state 3 C , 
antidiquark l dd in color state 3 C and noncolor state l ud , therefore the total color singlet can be 
constructed. 

The subamplitudes A\ (s, s 12345 , s 1234 , s l23 , s 12 ) and A\ u (s, Si 23A5 , s l23A , s 123 , s 12 ) are also the 
color singlets. 

The system of graphical equations Fig. 3 is determined using the selfconsistent method. The 
coefficients are determined by the permutation of quarks [44, 45]. We should discuss the coefficient 
multiplying of the diagrams in the equations of Fig. 3. For example, we consider the first subam- 
plitude A\ m (s , S12345, s i234, s i23, 512)- I n the Fig. 3 the first coefficient is equal to 2 (permutation 
particles 1 and 2). The second coefficient is equal to 6 = 2 (permutation particles 1 and 2) x3 
(we consider the third, the fifth, the sixth particles). The similar approach allows us to take into 
account the coefficients in all equations. 

In order to represent the subamplitudes A\ uu , A\ ud , A\ M , j^ uu ^ dd ^ ^r" u i Md i dd - m ^ e form of 
dispersion relations, it is necessary to define the amplitudes of qq and qq interactions. This is 
similar to the three quark case (Sec. II). 

We use the results of our relativistic quark model [38] and write down the pair quark amplitudes 
in the form: 

°n(«tt) = ; Gl !ff q , , (17) 

1 — tf n {Sik) 



(mj+m fc ) 2 A 



B n (s lk )= [ ^k P^ s 'ik)Gl{s' ik ) (lg) 

(mi+m k y 

Here G n (sik) are the diquark vertex functions (Table V). The vertex functions are determined 
by the contribution of the crossing channels. The vertex functions satisfy the Fierz relations. 
These vertex functions are generated from gy, gft and gy S \ B n (sik) and p„(sjfc) are the Chew- 
Mandelstam functions with cutoff A [34] and the phase spaces, respectively: 



/ jPC\ ( 1 jPC\ S ik . 01 jPC\ 1 r/ jPC\ ( m * m fc) 2 

pn{Sik,J ) = \a(n,J ) — | - + p(n, J ) + 5{n,J )- 



{mi + m k ^ s ik 



J (stk - (rrii + m k ) 2 )(s ik - (m 8 - m k ) 2 ) 
x . (19) 

Sik 

The coefficients a(n, J PC ), /3(n, J PC ) and 5(n, J PC ) are given in Table V. 

Here n — 1 coresponds to gg-pairs with J p = 0~, n = 2 corresponds to the gg-pairs with 
J p = 1~, n = 3 defines the gg-pairs with J p = + , n = 4 coresponds to J p = 1 + gg states. 

In the case in question the interacting quarks do not produce a bound state, therefore the 
integration in (20) - (24) is carried out from the threshold (rrii + rrik) 2 to the cutoff A(ik). 

The coupled integral equations correspond to Fig. 3 can be described as: 
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Al UU I \ \\B\uu{S\2) _ , -tUU\ A l uu I I \ 

^1 I s , S12345; s 1234, S123, S12J = 77 75 7 77 + 2Ji(Si2, 1 JA 1 (S, Sl2345j s 1234, ^123) s 13) 

[1 — #1««(S 12 JJ 

1 \ s i s 12345) s 1234) 5 123) 4) , (20) 

^1 IS, S12345, S1234, S123, S12J — 77 75 7 77 + ^H s 12, 1 J^l IS, S123455 s 1234, Sl23? S 13J 

[1 - 13^4512 J J 

+ 6Ji(Si2, l dd )A{ U \s, S12345, S1234, Sl23, «i 3 ) , (21) 



^2 ( S ) s 12345? S1234, S12, S34) 



^3 I s , S12345? s 12, S34, S56J 



Ai5 1U d(s 12 ) 

+ 2Ji(si2, r d )A} dd (s, S12345, S1234, S123, s'i 3 ) 



-Al («, Si2345,Sl234,Sl23,Sl 2 ) = 77 75 7 77 + 2Ji(Si2,l )^i (s, Si 23 45, Sl234, «123, S 13 J 

[1 — ±> ltt d(Sl2jJ 



+ 4Ji(si2, l ud )A\ ud (s, S12345, S1234, S123, s'13) 

+ 4J 2 (S12, r^Al^is, S 123 45, ^1234, s' 13 , S 24 ) , (22) 

A2-Bi«« (si2)-B 1 dd (S34) 
[1 - Si««(si 2 )][l - B l3d -(s u )] 

+ 2J 4 (S12, S34, 1"", l dd )A\ UU ( S , S12345, S1235, S125, S' 15 ) 

+ 2J 4 (si2, s 34 , 1 M , l uu )A\ dd (s, S12345, S1235, S125, sis) 
+ 4J 3 (si2, S34, l uu , l dd )A}" d (s, S12345, S1234, s' 123 , s 23 ) 

+ 4J 6 (si2, S34, 1 UU , r U )Al UUlM ( S , S12456, S1456, s' 15 , S«) 

+ 4J 8 ( Sl2 , S3 4, T", l^Ar 1 " 1 "^, *12345, S' 15 , 4, 4) , (23) 

A3-Bi«u (si 2 )-B 1U d (s 34 )i? 1 dd(s5 6 ) 
[1 - B 1 uu(s 12 )][l - S lud -(s 3 4)][l - B l3 d(s 56 )} 

+ 2J 9 (S 12 , S34, S 56 , I"", l^, 1 M )AI UU {S, S 123 45, S1234, si 23 , s' 2 



2:-! / 



+ 2Jg(s 12 , S34, S56, 1 , I' 5 Si2345> S1234, S 12 3, S23J 

+ 2Jg(s 12 , S34, S 56 , 1"", 1' 

j Si2345) S1234, S 12 3) S23) 

+ 4Jg(si2, S34, S56, 1"", 1' 

1 (s, Si2345) S1234, S 123 , S 23 ) 



where 



+ 2Jg(si2, S34, S56, 1" U )^4} (s, «12345, «1234, s' 123 , S 23 ) 

+ 4Ji (si2, S34, S56, 1"", l^)^" Uldd ( S , S12345, ^2345, 4 3 > S 4 5 ) , 

(24) 



(m 1 +m 2 ) Aj ^ 

w -x ^( g 12) / <fej 2 ^( S 12)P»( g/ 12) [ dZliX) , r\ 

Ji{s 12 ,i) = T —, TT / ; / — - — , (25) 

[1 - Bi(Si 2 )\ J n s' 12 - S12 /, 2 

(mi+m 2 ) 2 —1 



(m 1 +m 2 ) 2 A i 



■Msl2,0 



Gi(si 2 ) ? ds , 12 G i (s' 12 )p i (s' 12 ) 1 7^l(2) + /^2(2) 



[1 - BiM] J ir s' l2 - s 12 2tt 7 2 7 2 

(mi+m2j — 1 



-1 



^3(2)+ 

x / dz 3 (2) (26) 
23( V V 1 " - *Z(2) - z&2) + 2^(2)^(2)^(2) 



(m 1 +m 2 ) 2 A i 



f/ „ • ■ \ _ gi(si2)Gj(s34) / d S ' 12 G i (s' 12 )p l (s' 12 ) 



x 



[1 - Bi(s 12 )][l - Bj(su)] , J,* s' 12 -s 12 

(m.i+m.2) 

(m 3 +m 4 ) 2 Aj 

rf4 4 G , i (s / 34 )p i (s / 34 ) /■ d^i(3) /■ dz 2 (3) 



7 7T s'o 4 -S34 7 2 7 2 ' 1 ' 



s 34 — s 34 

(m3+ni4) 2 —1 —1 



^4(312,334,1, j) = ^ ^^^)j 7l(gl2,0 , (28) 
7 6 (S12,S34,«, j) = Jl(si2,i) • 7l(s 3 4, j) , (29) 



7 8 (Sl2,S34,^ j) 



(m 1 +m 2 ) Aj 

GiWGiW ~~ f ds' 12 G t (s' 12 ) Pl (s' 12 ) 



[1 - Bi(s l2 )][l - Bj(s 3A )] so ^* S12-S12 

(m.l+m.2) 2 



(m 3 +m4) 2 Aj 



X 



X 



X 



/ * 



C?S 3 4 ^( s 34)Pi( s 34) 



(m 3 +m4) 2 34 

j_ 7^(8) 7^(8) 7^(8) Z4 f\ 7^(8) 26(8)H 

(27T) 



-1 -1 -1 z 4 (8)- -1 
1 

1 - ^(8) - ^ 2 (8) - zf(8) + 2^(8)^3(8)^4(8) 



x 1 (30) 

- ^2(8) - 4(8) - zj(8) + 2z 2 (8)z 5 (8)z 6 (8) 
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^9(S12,S34, S56,2, j, k) = j^fgL J 3 (s 12 ,su,i,j) , (31) 

[1 - t>k{Sb§)\ 



J\o{si2i S34, s 56 , i, j, k) — 



G ! i (si2)Gj(s34)Gfc(s5 6 ) ? <fe' 12 G , i (s' 12 )pj(s / 12 ) 



[1 - ^( Sl2 )][l - ^( S 34)][l - 5 fc ( S56 )] J TT S' 12 - S 12 



(m 3 +m 4 ) 2 A J - (m 5 +mg) 2 A fc 



J TT S04 — S34 7 7T 



+ 1 +1 +1 +1 25(10)+ 

1 / dzi(10) /■ ^2(10) /■ dz 3 (10) /■ ds 4 (10) 



x 



-1 -1 -1 -1 *5(1-)- 

x 1 (32) 

^/l - z?(10) - zf (10) - zl(10) + 2zi(10)z 4 (10)z 5 (10) 

In the equation (25) Zi(l) is the cosine of the angle between the relative momentum of particles 

1 and 2 in the intermediate state and the momentum of the particle 3 in the final state taken in 
the cm. of particles 1 and 2. We can go from the integration of the cosine of the angle dzi(l) to 
the integration over the subenergy ds' l3 . 

In Eq. (26) 21 (2) is the cosine of the angle between the relative momentum of particles 1 and 

2 in the intermediate state and the momentum of the particle 3 in the final state taken in the 
cm. of particles 1 and 2, -22(2) is the cosine of the angle between the momenta of particles 3 and 
4 in the final state of cm. of particles 1 and 2, 23(2) is cosine of the angle between the relative 
momentum of particles 1 and 2 in the intermediate state and the momentum of the particle 4 in 
the final state of cm. of particles 1 and 2. Then we pass from dzi(2)dz2(2)dzs(2) to Gk' 13 <is3 4 <is 24 . 

In Eq. (27) 21 (3) is the cosine of the angle between the relative momentum of particles 1, 2 
in the intermediate state and the relative momentum of particles 3, 4 in the intermediate state 
in cm. of particles 3 and 4; ^(3) is the cosine of the angle between momentum of particle 3 in 
the intermediate state and relative momentum of particles 1, 2 in the intermediate state in cm. 
1 and 2. We pass from dz 1 (3)dz 2 (3) to ds'^ds^. The similar method is used for the functions 
(28), (29), (31). 

In Eq. (30) 24(8) is the cosine of the angle between momentum of particle 5 in the final 
state and the relative momentum of particles 1, 2 in the intermediate state in cm. of particles 
1 and 2; 2:2(8) is the cosine of the angle between the relative momentum of particles 1, 2 in 
the intermediate state and the relative momentum of particles 3, 4 in the intermediate state in 
cm. of particles 3 and 4; 2:3(8) is the cosine of the angle between momentum of particle 3 in 
the intermediate state and the momentum of particle 5 in the final state in cm. of particles 1 
and 2; 24(8) is the cosine of the angle between the momentum of particle 3 in the intermediate 
state and the relative momentum of particles 1, 2 in the intermediate state in cm. of particles 
1 and 2; 2 5 (8) is the cosine of angle between momentum of particle 6 in the final state and the 
relative momentum of particles 1, 2 in the intermediate state in cm. of particles 3 and 4; 2 6 (8) 
is the cosine of the angle between momentum of particle 6 in the final state and the relative 
momentum of particles 3, 4 in the intermediate state in cm. of particles 3 and 4. We pass from 
dz\ (8)^22(8)^23(8)^24(8)^25(8)^26(8) to ds'^ds'^ds'^ds^ds'^eds^Q. 

In Eq. (32) 21 (10) is the cosine of angle between relative momentum of particles 1, 2 in the 
intermediate state and the relative momentum of particles 3, 4 in the intermediate state in cm. 
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of particles 3 and 4; 2:2 (10) is the cosine of angle between the relative momentum of particles 1, 
2 in the intermediate state and momentum of particle 3 in the final state in cm. of particles 1 
and 2; 23(10) is the cosine of the angle between the relative momentum of the particles 3, 4 in the 
intermediate state and the relative momentum of particles 5, 6 in the intermediate state in cm. 
of particles 5 and 6; 24(10) is the cosine of angle between relative momentum of particles 1, 2 in 
the intermediate state and the momentum of particle 5 in the final state in cm. of particles 3 
and 4; 25(10) is the cosine of the angle between the relative momentum of the particles 3, 4 in the 
intermediate state and the momentum of particle 5 in the final state in cm. of particles 3 and 4. 
We pass from ^21(10)^22(10)^23(10)^24(10)^25(10) to ds'^ds^ds'^ds'^ds'^. 

Let us extract two- and three-particle singularities in the amplitudes A\ uu (s, S12345, S1234, S123, S12), 

jl ^dd / \ .^ud, \ .-^uu^dd , \ 

A \ [S, S12345, S1234, S123, S 12), A \ \ s i s 12345, S1234, S123, S 12j) A 2 (5,512345,51234,512,534), 

.^uu-^ud-^dd / % 

^3 (5,512345,512,534,556): 



(33) 



[1 


— B]y.u (s 12 )] 




1 dd / 

«i (s, S12345, 


51234, Sl23, Sl2)-B 1 dc 


K S 12) 


[1 






1 rid , 

«1 (5, S12345, 


51234, Sl23, Sl2)-B lu , 





A \dd , , OL x (S, S12345, 51234, 5123, 5l2j-D 1 dd(Si2j , 
A \ (5, S12345, 51234 , 5l23, 5l2j — f. 75 7 yj , (<j4J 



^1 (5,512345,51234,5123,512) — 7\ 5 7 yj , [OO) 

[i — n^dySu)] 

iuui dd. . a\ 1 (s, S12345, 51234, S12, S34)5i«u(s 1 2)-B 1 ,id (S34) , . 

A 2 [S, S12345, 51234 , 5l2, S34J - r- — , (dOj 

[1 - Si^(s 12 )J[l - B lSS {s M )\ 

luulu d 1 dd, . Oi\ 1 1 (s, S12345, 5l2, 534, S56)-Bl««(si2)-B lu d(s34)i? 1 dd(s56) 



^3 (5, S12345, 5i2, 534, S56J 



[1 - B X uu (s 12 )] [1 - B lu d(s 34 )] [1 - ^dd (5 56 )] 

(37) 



We do not extract four-particles singularities, because they are weaker than two- and three- 
particle singularities. 

We used the classification of singularities, which was proposed in paper [39]. The construction of 
the approximate solution of Eqs. (33) - (37) is based on the extraction of the leading singularities 
of the amplitudes. The main singularities in Sik = (m ; + m k ) 2 are from pair rescattering of the 
particles i and k. First of all there are threshold square-root singularities. Also possible are pole 
singularities which correspond to the bound states. The diagrams of Fig. 3 apart from two-particle 
singularities have triangular singularities and the singularities defining the interactions of four, five 
and six particles. Such classification allows us to search the corresponding solution of Eqs. (20) 
- (24) by taking into account some definite number of leading singularities and neglecting all the 
weaker ones. We consider the approximation which defines two-particle, triangle and four-, five- 
and six-particle singularities. The contribution of two-particle and triangle singularities are more 
important, but we must take into account also the other singularities. 

The five functions «j are the smooth functions of s ik , s^-fe, s^/m as compared with the 
singular part of the amplitudes, hence they can be expanded in a series in the singularity point 
and only the first term of this series should be employed further. Using this classification, one 
defines the reduced amplitudes as well as the _B-functions in the middle point of physical region 
of Dalitz-plot at the point Sq: 
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s + 4 £ mf 

so = ~^r^- , (38) 
E < 

i<fc 

3 3 

S123 = S E m ifc - E ( 39 ) 
»,fc=i i=l 

i<fe 



Sl234 = S E m 'fc- 2 E m ?- ( 40 ) 

i,fc=l j=l 

Such choice of point so allows us to replace integral equations (20) - (24) (Fig. 3) by the 
algebraic equations (41) - (45), respectively: 

«r = A + 2/ 1 (rr M )af + 6/ 1 (ri ,Mi ")af , (41) 

af = X + 2I 1 (l d h dd )af + 6h(l d h ud )al u \ (42) 

= A + 2/ 1 (rt u >r + 27 1 (l ud "l^> 1 ld " d " + 4/ 1 (l ud > d ~)< <? + 4/ 2 (r d >n^>r id ~ J , 

(43) 



= A + 2I A {l uu l M r u )ar + 2h{l d h uu l dd )af + Ah^l^r^a 1 ^ 



+ Ah{i uu i dd i uu i dd W; ui + Ah {r u i dd r u i ud \ dd )ar i 1 , 

(44) 

= A + 2 / 9 (i»f' i "l JJ f >f + 2/ 9 (l (i " J f J l»l d "Vl S + 2/ 9 (i Mu r J i d "' i r J )a! u ' 



+ 4/ 9 (r M i JJ r J r J )a^ d + 2/ 9 (i JJ r J r M i Md ")a^ d + 4i 10 (r u rh dd r u i dd )al uuldd , 



(45) 



where Aj are the current constants. We used the functions h, I 2 , I3, h, h, h, hi ho'- 



B^) f ds' 12 Gj{sf) Pl {s' 12 ) + rd Zl {l) 1 



m) B^stf) I tt \' l2 -sf J 2 l-S^)' (46) 

(mi+m2) — i 



J 2 (yfc) 



B s ( 8 ^)B k (^) f ds' 12 G 2 M 2 )p t (s' l2 ) 1 + /^(2) + /^ 2 (2) 



^(sq 2 ) , J* s' 12 -sf 2nJ 2 J 2 



(mi+m2) 2 — 1 — 1 

^3(2) 

X 

*s(2) 



/ <** 3 (2) 

.(V V 1 - ^) - z 2 2 (2) - 4(2) + 2^(2)^(2)^(2) 



12 



h(ijk) 



I 8 (ijklm) 



I w (ijklm) 



l- Bj {s^)l-B k {s' 2i )' 

(m 1 +m 2 ) 2 A i 

B k (sf) f ds' 12 G}{sf) Pi {s' 12 ) 



(m3+m 4 ) 2 A J - 



X 



ds' zi G]{sf) Pj {s'^) 7^(3) + r 1 dz 2 (3) 



s' u - sjj 4 J 2 J 2 1 - B k (s' 23 ) ' 

(m 3 +m4) 2 "I -1 



/" Qggi^jifo J Pi { S 34> f 

J 7T — J 

+ 1714 

h{ijk) = h(ik), 
I 6 (ijkl) = h(ik) ■ h{jl) 



(m 1 +m 2 ) 2 At 



B^B^BUsf) f ds> 12 GlisfWu) 



x 



Bi{sf)B 3 {sf) J tt S ' 12 -4 2 

(mi+m2) 

(m 3 +m 4 ) 2 A J - 

~ ds' u G]( s ^) Pj ( s ' u ) 



J 71 



S34 Sg^ 
(m3+m,4) 2 



+ 1 +1 +1 24(8)+ +1 z«(8) + 

1 + f d Zl (S) + f dz 2 (S) + f dz 3 (&) f ) <fe 5 (8) f 

X WJ—J—J— J dz ^J— J iz ° 

-1 -1 -1 24(8)- -1 26(8)" 



1 



X 



X 



X 



1 - ^(8) - zl(8) - zj(8) + 2^(8)^(8)^(8) 
1 

y/l - 5(8) - 5(8) - ^ 6 2 (8) + 2z 2 (8)^ 5 (8)^ 6 (8) 
1 1 1 



1 - S fc (s' 15 ) 1 - 5,(4 3 ) 1 - S m (4 6 ) ' 
ig(u'A;Z) = h(ijl) , 



(m 1 +m 2 ) 2 A i 



B^BUsf) f ds' 12 GZ(sl 2 ) Pi (s> 12 ) 



x 



B^Bj^B^) J tt S ' 12 -^ 2 

(mi+m 2 ) 2 

(m 3 +m 4 ) 2 Aj (m 5 +mg) 2 A fc 

~~ 7 ^34 ^(^)^(44) > ^ 6 ^(^ 6 )p fe (4 6 ) 

J 71 S34 Sg J 7T Sgg Sg 

(m 3 +m 4 ) 2 (m 5 +m 6 ) 2 
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X 




-1 -1 -1 -1 *5(1-)- 



1 



X 



yji - zfiio) - 4(io) - 5(io) + 2^(10)^4(10)^5(10) 
1 1 



X 



1-^(43)1-^(545) ' 



(53) 



where i, j, k, I, m correspond to the diquarks with the spin-parity J p = + , 1 + and mesons with 
the spin-parity J p = 0~, 1~ . 

The other choices of point s do not change essentially the contributions of ctj, therefore we 
omit the indices sff. Since the vertex functions depend only slightly on energy, it is possible to 
treat them as constants in our approximation. 

The solutions of the system of equations are considered as: 



where zeros of D(s) determinants define the masses of bound states of baryonia. 

As example, we consider the equations for the quark content uuuddd with the isospin 1 = 3 
and the spin-parity J p = 3~ (Fig. 3). The similar equations have been calculated for the isospin 
I = 0, 1, 2, 3 and the spin-parity J p = - , 1~, 2~, 3~. We take into account the u and d quarks. 

In Appendix I the reduced amplitudes of baryonium uuduud 1.1 = 00 are given. 

IV. Calculation results. 

The poles of the five reduced amplitudes correspond to the bound state and determine the 
mass of the baryonium with the quark content uuuddd, with the isospin 1 = 3 and the spin-parity 
J p = 3-. The quark mass of model m = 410 MeV coincides with the ordinary baryon one in our 
model (Sec. II). 

The model in question has only two parameters: the cutoff parameter A = 11 (similar to the 
three quark model (Sec. II)) and the gluon coupling constant g = 0.314. This parameter is 
determined by the baryonium mass M = 1835 MeV . The estimation of theoretical error on the 
baryonia masses is 1 MeV . This result was obtained by the choice of model parameters. 

We predict the degeneracy of the some states. In the Table II the calculated masses of non- 
strange baryonia are shown. The contributions of subamplitudes to the six-quark amplitude are 
shown in the Appendix I (for example, the baryonium with the mass M = 1835 MeV). The states 
(AA + An + nA + nn) and (AA + Ap + pA + pp) with the isospin 1 = and the spin-parity 
J p = 0~ possess the mass M = 1835 MeV. For the (AA + nA) and (AA + An) with the isospin 
1=1 and the spin-parity J p = 0~ we obtained the mass M = 1928 MeV. 

We predict the degeneracy of baryonia M(uudddd, I = 2) = M(uuuudd, I = 2) = M(uddddd, I = 
1) = M(uuuuud, I = 1). For the states M(uududd, I = 1) = M(uddudd, I = 0) = M(uuduud, I = 
0) and M(uuuuuu, I = 0) = M(dddddd, 1 = 0) the degeneracy is also obtained. 

V. Conclusion. 

A somewhat simple picture of baryonium is that of a deuteron-like NN bound state or reso- 
nance, benefiting from the attractive potential mediated by the exchange of gluon [33] . We do not 
consider the influence of annihilation on the spectrum. 

Entem and Fernandez, describing scattering data and mass shifts of pp levels in a constituent 
quark model, assign the threshold enhancement to final-state interaction [46, 47]. Zou and Chi- 



on(s) = F J (s,A l )/ J D(s) 



(54) 
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ang find that final state interaction makes an important contribution to the pp near threshold 
enhancement [48]. 

The baryonium state with M = 1835 MeV is considered as pp state [3] or the second radial 
excitation of rf meson [49]. 

In our case this state have following content AA + Ap + pA + pp with isospin 1 = and 
spin-parity J p = 0~ . 

We calculated the masses of baryonia with the isospin / = 0, 1, 2, 3 and spin-parity J p = 
0", 1", 2", 3" (Table II). 

The interesting reseach is the contribution of baryonia consisting of u, d, s-quarks and anti- 
quarks. 
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Appendix I. The reduced amplitudes of baryonium IJ = 00(1835). 



af = A + 2aOi(l u "l u ") + 2 af h(l uu l ud ) + 2 af h(l uu ud ) + i af h(l uu uil ) 

+ 2a1 ud I 1 (l uu ucl ) 

af = A + 2cOi(l""l u ") +2of/ 1 (ri*) + 2af d /i(l fifi fiJ ) + 4 af ^(l^O"") 

+ 2af/i(l 6S ,,B ) 

a}" = A + afAll™!™) + «r^i(l Mfi l™) + 2ar^i(l ufi l ufi ) + afh(r u r d ) + af h(r u l m ) 

+ afh{l u % ud ) + af d h{r^ d ) + 2 af ° ^(rV) + af ^(r^) + af h{r^) 

+ a^/^rro^) + a f " d ; 2 (i-r s o U(i ) + af ^^(r^o^) 

af = A + cOi(rt u ") + «r S A(l u<i "l B ") + af^l"^) + afh(rh dd ) + af h(l ud ~0 ud ) 

+ 2 af ^(rt^) + 2 af * h(l ud Q uU ) + af ^(rt^) + af'/^l^O^") 

+ 2 ar^/^l^l-O^) + 2 af ^(l^O^O**") 

af = A + af"/i(l d "l™) + arA(l d "l u ") + af Ji(l ds l d ") + af h(l m l dd ) + 2 af h{l dn $ ud ) 

+ af d I 1 {l M Q ud ) + 2 af °/i(l <ia uO ) + af h(l du du ) + af J /i(l ds <w ") 

+ 2 a^° ud / 2 (l ds 1^0" d ) + 2 af J 2 (l da ud fid ") 

a}" = A + afh(l dd r d ) + afh(l dd l m ) + 2 afh(l dd ud ) + 2 af h(l dd ud ~) 

+ 2 afh(l dd ud ) + 2 afh(l dd ~0 m ) + 2 af^h (l<%"*0^) 

af = A + af"/ 1 (0 ud l uu ) + 2 af h(0 ud V u ) + af /^O^) + 2 afh(0 ud l m ) + af/^O^) 

+ af'A (0 ud ud ) + 2 af "/i(0 ud u ") + af *h (0 ud uJ ) + 2 af /^O^") + af h(0 ud dd ) 

af = A + af/i (0" J 1™) + 2 af 7i ((f d >") + 2 af 7i (0 sJ l uJ ) + af h (0^1 ds ) + af 7i (0" d "l dJ ) 

+ af A (0" J 0" J ) + 2 af/i (0^0™) + 2 af h (0"V d ") + af ^ (0 fiJ dfi ) + af A (O^O**) 
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af ° = A + ai"Vi(O uS l uu ) + afh(O ua l m ) + 2 af h(O uU l uU ) + af 7 1 (0 ufi l ud ") + af ^((T^) 

+ af d 7i(0 u "0 ud ) + af d "7i(0 u "0" J ) + 2 af" Ji(O ua O ua ) + af d ~7i(0 ua uJ ) + af 7 1 (0 ua da ) 

+ ^""^^(OTr) + af°* d ~7 2 (0 ua r u a<i ~) + ar° ud 7 2 (0 M "l""0" d ) + «f d0 " d 7 2 (0 ua nd ad ~) 

a?" d = A + af "A (0 ud > u ) + 2 afA (0 ud > a ) + afh (0 ud l ud ) + afh (O^l**) + afh (0 ud ud ) 

+ 2 afh (0 uJ aJ ) + 2 af^ (0 uJ ua ) + af ^ (0 ud uJ ) + afh (0 ud dd ) 

+ af^h^h^ 3 ) + 2 af "7 2 (0 uJ U(i aJ ) 

af = A + af 7 1 (0 da l aa ) + 2afh(0 du r u ) + af h(O du l du ) + afh(O du l dd ) + 2 af h(0 du ud ) 

+ afh(0 du ud ) + 2 af h(0 du uii ) + af 7 1 (0 da da ) + afh(0 du dd ) + af °" d 7 2 (0 da l aa ud ) 

+ 2af d ° M > 2 (0 da w a < 1 ) 

af = A + 2 afh(O dd lf + 2 afh{^l m ) + 2 a?" ^(0^) + 2 af ^ (0 dd ~0 ad ~) 

+ 2afh(0 dd ud ) + 2afh(0 dd ~0 du ) + 4 af°" I 2 (0 dd ud ud ) 

»„!«. = A + 2 ^j^^Sq^ + 2 a f J^!™^™^) + 4 a 0^J 3 ( 1 «« 1 «« ™) 

+ 4 a° ud ° M 7 6 (l M "l""0" d 0" J ) + 4 ajf ^ ffl ^J 8 (l Utt l aa 0^0 ua 0**) 

««o«* = A + a i«« j 4 ( **i««i™) + 2 af ^(i^o^'r") + 2 af 7 3 (l" M a<7 l Md ~) + 2 af d h(l uu O nd O ud ) 

+ af "/ 4 (O fiJ l n "0" J ) + 2 af ^7 3 (l" M 0" (! n ") + 2 af d 7 3 (l" M ad ~0 Md ~) + 2 af° ud I 6 {r u ud ud l m ) 

+ 2 af° ud l 6 (l uu ud ud ud ) + 2 Q ,i^i" d > d / 8 (i««o" J O ud l ud "r") 

+ 2 af "O^J^UUQUdQUdQUUQUd^ 

u *° ud = A + af 7 4 (O ud l aa l uu ) + 2 af 7 3 (l aa Mc( l Ma ) + 2 af 7 3 (l aa Md l da ) + af h(0 ud l m ud ) 

+ 2 af h{l m ud ud ) + 2 af I 3 (l ua ud uU ) + 2 af l 3 (l m ud du ) 

+ 2 ar°"l 6 (l uu ud ud r u ) + 2 af ^(fWo"") + 2 a^^l^O^l* !™) 
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a f d ° ud = A + arh(O ud O a h uu ) + af* h(O ud O ud l m ) + af* I 3 (0 ud a h uU ) + af I 3 (0 ud 0^1 ud ) 

+ a{ m l 3 (0 ud ad l du ) + af l 3 (0 ud ad l dd ) + af d h(O ud O ud O ud ) + af ^ h(O ad O ud O ud ) 

+ af* I 3 (O ud O ud O uii ) + af " I 3 (0 ud ad ud ) + af u I 3 (0 ud ud du ) + af ' I 3 (0 ud ad dd ) 

+ a 1 ™ 1 ** 7 6 (0 ud fid ~l uu l fifi ) + af u ^ , J 6 (O^O i *r u O fid ~) + af * °" d / 6 (0 u,i fiJ u,i l fifi ) 

+ af d °^ 6 (0 ud sd > d (F J ) + af "^"^(O^rV^) + h{^ ad l uu ^ ad ) 



al™!™** = A + 2 af d / 9 (l uu <w ~l fifi l ud ~) + 2 af fl / 9 (l fflS d,i "l uu l ds ) + 2 af d / 9 (l uu <w ~l fifi (r i ) 

+ 2 af* J 9 (l fflffl < "l uu aJ ) + 4 af V 9 (l uu l""0 <iJ O u ") + 2 af d / 9 (l UM dd ~r"0" J ) 

+ 2 a°f l 9 (l m dd l uu du ) + 4 af" ^ /^(l^O^l^O^O^) 

ai^i^o^" = A + a i <H, / 9 (i*»o fiJ l uw l iia ) + 2 af V 9 (l Wtt fid ~l <ia l uii ) + 2 af J J 9 (l ttU fiJ l <ia l ,ld ~) 

+ af* I 9 (l du O ud r u l m ) + a{ di I 9 (l du ad r u l dd ) + 2 af" l 9 (l uu l du ad ud ) 

+ af d "/ 9 (i^o" d "l uu O" J ) + af * (2 / 9 (r«l d «0" J u ") + 2 / 9 (i««o^l d "0™)) 

+ 2 af ^(i^o^l^O"**) + af* I 9 (l du O ud l uu O m ) + af l 9 (l du ud l uu dd ) 

+ 2 af fl °" d 7 10 (i««i^ " J u<i r") + 2 af dolid / 1 o(l M "l d "O sJ 0^0" J ) 

+ 2 af o / 9 (l fifi ud l uJ l dfi ) + af I 9 (l ud ud l aa l dd ) + af d / 9 (l" J 0" d l™0 Mc( ) 

+ 2 af l 9 {l m l ud ud ud ) + af °(2 ^(l^l^O^O" ) + 2 / 9 (i«« MC 4 M<7 M ")) 

+ af d / 9 (l M<7 Md l""0" J ) + 2 af %(l""O ud l u<i ~O d ") + af d / 9 (l uJ ud l fifi dJ ) 

+ 2 af "^"/^(i^i^o^O^r") + 2 af ^^(l^l^rt^O^) 

af d °" " d = A + a{ uu I 9 (0 ud uU ud l uu ) + af I 9 (0 ud uU ud l uu ) + a}" (I 9 (0 ud uU ad l uU ) 
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a 1 ^ (I 9 (0 ud uU ud l du ) + I 9 (0 ud ad uU l du )) + af I 9 (0 ud ud uU l dd ) + aT" h(O ud O ua O ad O ud ) 
af d h(O ad O uU O ud O ud ) + aT* (l9(0 ud uU ad uU ) + I 9 (0 ud uU ud uU ) + I 9 (0 ud ud uU uU )) 
a?" \l 9 {0 ad uU ud ud ) + I 9 {0 ud ad uU ud )) + a ^ {I 9 {0 ud uU ud du ) + I 9 {0 ud ad uU du )) 
af d I 9 (0 ud ud uU dd ) + al uuiau I w (O ud O uU O ud l uu l uu ) + a^ ^ I 10 (O ud O uU O ud l uu O ud ) 

^ J^QUdQUUQUdQUd^ + a O»^ /io ( «d «« «d «d tid) 
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Table I. S-wave baryon masses M(J P ) (GeV). 





V 2 / 






V 2 / 




N 


0.940 
0.940 


(0.940) 


A 


1.232 
1.232 


(1.232) 


A 


1.022 

1 .U»o 


(1.116) 


£* 


1.377 

L.Ol I 


(1.385) 


£ 


1.050 
1.193 


(1.193) 


"* 


1.524 
1.524 


(1.530) 




1.162 
1.325 


(1.315) 




1.672 
1.672 


(1.672) 



Table II. S-wave baryonia masses. Parameters of model: cutoff A = 11.0, gluon coupling 
constant g = 0.314. Quark masses m U: d = AlQMeV. 



I 


Quark content (baryonia) 


J 


Mass (MeV) 





uuuuuu (AA), 
dddddd (AA) 



1 
2 
3 


1973 
1824 
1938 
2290 




uddudd (AA + An + nA + nn), 
uuduud (A A + Ap + pA + pp) 



1 
2 
3 


1835 
1784 
1851 
2455 


1 


uddddd (AA + nA), 
uuuuud (AA + Ap) 



1 
2 
3 


1928 
1770 
1857 
2395 




uududd ( AA + An + pA + pn) 



1 
2 
3 


1835 
1784 
1851 
2455 


2 


uudddd (AA +pA), 
uuuudd (AA + An) 



1 
2 
3 


1928 
1770 
1857 
2395 


3 


uuuddd (AA) 



1 
2 
3 


2067 
1783 
1938 
2290 
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Table III. 13 = 00, uuduud (1835 MeV), A = 11.0, g = 0.314. 



Snbamnli tildes 


Contributions nerrent 




3.7 




3.7 




9.0 


j^ud 


7.8 




7.8 


j^ydd 


6.7 


^QUd 


3.6 


j^ud 


3.6 




7.0 


^QUd 


6.8 




6.8 


^Qdd 


6.6 


A 1 UU 1 uu 

A\ 1 _ 


2.4 


. -^uuQud 

A 2 


2.0 


A 2 


2.6 


, ud ud 

^-2 


2.8 


a uu ^dd uu 

A 3 


4.0 


a yuu -^duQud 

A 3 


4.5 


, -^uu -^udQud 
^3 


4.5 


^QudQUU Qud 


4.3 


EA 


73.0 


E^2 


9.8 


E^3 


17.2 



Table IV. I J = 33, uuuddd (2290 MeV), A = 11.0, g = 0.314. 



Subamplitudes 


Contributions, percent 




9.9 


j^dd 


9.9 


j^ud 


25.4 


jl ^uu^dd 

A 2 


14.5 


. -yuu -ydd -yud 

A 3 


40.3 


E^i 


45.2 


E^2 


14.5 


E^3 


40.3 
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Table V. Vertex functions and Chew-Mandelstam coefficients. 



i 


{Ski) 




R- 
Pi 




+ diquark 
f + diquark 
0~ meson 
1~ meson 


Ag Sgm 2 kl 
3 (3s fci ) 

2<? 


1 
2 

1 


1 (m fc -r» ; ) 2 

2 (m k +mi) 2 

4m k mi 1 




1 (m k -mi) 2 


3 

3 (3s fc; ) 
43 


3 

1 
2 

1 


3(m fe +m i ) 2 6 

1 (m fc -m ; ) 2 

2 K+m,) 2 

im k mi 1 


6 (m fc +m ; ) 2 



1 (m k ~rni) 2 


3 


3 


3(m fe +mi) 2 6 


6 (m k +mi) 2 
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+ 




Fig. 1. Diagrams which correspond to a) production of three quarks, b - f) subsequent pair 
interaction. 




25 




26 




27 




Fig. 3. Graphic representation of the equations for the six-quark subamplitudes A\ (I = 1,2,3) 
in the case of baryonium uuuddd I J = 33. 
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